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ABSTRACT 


The instability of the ordinary electromagnetic mode propagat- 
ing perpendicular to an external magnetic field is studied for a 
single-species plasma with ring velocity distribution. The marginal 
instability boundaries for both the purely growing mode and the 
propagating growing modes are calculated from the instability cri- 
teria. The dispersion characteristics for various sets of plasma 
parameters are also given. The typical growth rates are of the 
order of the cyclotron frequency and enhanced by increasing P . 
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I. INTRODUCTION 

The instability of the ordinary electromagnetic mode propagat- 
ing perpendicular to an external magnetic field has received consid- 
erable attention. Hamasaki ( 1968 a, b) and Davidson and Wu (1970) 
investigated a stationary plasma with anisotropic temperatures and 
showed that for temperature anisotropy T /T ^ < 1 and sufficiently 
large p > 2 , where p |( is the ratio of parallel kinetic pressure to 
magnetic pressure, purely growing ordinary electromagnetic modes 
exist for propagation perpendicular to the external magnetic field. 
Counter-streaming plasmas with and without temperature anisotropy 
have also been previously investigated by many authors for the ordi- 
nary mode instability (Gaffey, Thompson and Liu, 197 2 , 1973; Lee 
and Armstrong, 1971; Bornatici and Lee, 1970; Tzoar and Yang, 1970). 

In this paper we wish to present a comprehensive picture of the ordi- 
nary electromagnetic mode for a single-species plasma. From the 
study of the dispersion relation derived from the linearized Vlasov- 
Maxwell equations, it is found that the instability can occur for 

6 > 1 and is dependent on the values of P and a, where P is the 
ll J- J- 

ratio of perpendicular kinetic pressure to magnetic pressure and a 

p 

is the ratio v /c , with c as the velocity of light, v is the aver- 

Oj/ 01 

age perpendicular particle velocity. Attention is called to the pre- 
diction of propagating instability with the real part of the complex 
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frequency, <u r , being non-zero. A detailed analysis is made of the 
steady state conditions of the plasma to determine the stability 
boundary. The solution of the very complicated dispersion relation 
is done numerically and compared with the analytical solutions that 
are obtained in certain limiting cases. The justification of the 
instability boundaries predicted by the instability criteria which 
require the instability to be absolute is also made through the 
numerical solutions of the dispersion relation. The growth rates are 
found to be typically of the order of the cyclotron frequency. 

The electromagnetic instability is important in high-(3 plasmas 
and is of interest in a broad range of astrophysical and laboratory 
applications. Recently the experimental observation of the electric 
field fluctuations made by the IMP-6 spacecraft in the outermagneto- 
sphere was reported by Gurnett and Shaw (1975). They observed strong 
electromagnetic radiation of harmonically related bands in the 
frequency spectrum with electric field fluctuation parallel to the 
external magnetic field. They suggested that it is due to an ordi- 
nary electromagnetic plasma instability and this suggestion is just- 
ified by our work. ■ 

The basic plasma model that is investigated consists of an 
equal number of electrons and ions immersed in a uniform external 
magnetic field and free of any electric field. Ions are considered 
to form a uniform positive background. In section II the dispersion 
relation is obtained. Various limiting forms of the dispersion 
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relation are discussed and general stability considerations are given. 
Section III presents the analysis of instability. Three conditions 
for instability are given and the instability boundaries are calcu- 
lated. In Section IV the numerical solutions of the dispersion 
relation of the various parametric dependences are discussed. 

Finally the conclusions of this study are summarized in section V. 
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II. DISPERSION RELATION 

For an infinite collisionless single-species plasma immersed 
in a uniform, static magnetic field ? we take the equilibrium dis- 
tribution function to be an even function of the component of 

velocity parallel to so that the ordinary mode decouples from 

^ 

the extraordinary modes for wave vector k j. B (k = k£ and 

O X 

B s B S )■ The ordinary mode is a purely electromagnetic wave 
o o z' 

with the electric field E |[ B^. 

The electric field in the (&,<ju) space for the ordinary mode 
is (e.g., Baldwin, Bernstein and Weenink, 1969) 

*<>-> ■ (1) 


where the numerical function is 


N(k,ti)) 


2 

(IXI) r 

= itue + ic£ x 3 + & z j dv v^ g(k, v^, v ( ) 


A J(p) 

i - J „ (uj + nm ) 
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2 . 
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- 6 -2— 

z me 


f / &f &f \ a nJ ( 

I 0 0 sp n 

J av v « r» ar - T i s^JJL ~ r 


sfj s. g; J n (p ) 

cu + nu> c ) 


( 2 ) 
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and the dispersion function is 


D(k,tu) = uj c 


2 , 2 2 

; k - u) + 


“ A n (k)ntu c 


n=-«> 


a) - n<n 


( 3 ) 


A„(fc) = 2 ” w p J n^P^ v n 


n 


5f 

dv dv 

Sv || i 


p p 

with nT = l4rrn e /m; m = eB /me; p = kv /cu , ; J is the nth order 

O C O 1 v U 

Bessel function of the first kind; f is the equilibrium distribu- 


tion function; and e, b, g are defined by 


(e, t, g) = f-^-x e"^* X [E (t = 0), B (t = 0), f (l) (t = 0) ] . 
J (2n Y 

The singularities at <u ~ nuo c in N(k,cu) are matched by the 
sing ular ities at u) = nou c in D(k,u>). Therefore all of the singularities 
in E(k,uo) come from the zeros of the dispersion relation D(k,ui) = 0. 

If we solve for u> (k real), exponential growth and dancing are implied 
for > 0 and uk < 0, respectively, while is the imaginary part 
of (1). 

In general the infinite sum in Eq. (5) converges rapidly for 
arguments of the Bessel functions less than order unity because 
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J n (p) “IFjrlf)' ' for p «1 . 

However for p larger than order unity the series converges very 
slowly because 


“ jj 0 ° s ( p -?-?) o » 1 • (5) 

i 

Therefore a large number of terms must be kept in order to obtain 
accurate results. This extensive calculation can be avoided by 
using both the integral representation of the Bessel functions and 
the Sommerfeld-Watson transformation technique (see Appendix A). 
The dispersion function then reduces to 


D(k,<i>) 


2 

GO 







sin fin 



dr sin fiT sin t 


J 

o 



T 

cos 2 


(6) 


where fi = o)/^. Clearly this expression is defined for all fi except 
possibly at ft = n and hence is the proper analytic continuation of 
Eq. (3). This form is particularly convenient for computational 
purposes. Before solving the dispersion relation exactly, it is 
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important to examine some interesting features of the dispersion 
relation. In the limit of a cold plasma system, the dispersion rela- 
tion becomes D(k,<jb) = <d 2 - c 2 ^ 2 - u> 2 = 0, which shows how the 
presence of the plasma alters the vacuum propagation of an electro- 
magnetic wave. In particular, u> is independent of B q . In the 
cutoff limit (k -> 0), the dispersion relation has roots at the 
plasma frequency and all harmonics of the cyclotron frequency except 
n = 0. In the resonance limit (k -+ «), the dispersion relation has 

roots at (y = ± (c^k. 2 + and uo = n<u c (n = ±1, ±2, ...). In 

the l imi t of very strong B , <d -* 00 and hence p -» 0, the dispersion 

relation has roots at k = ± (c^ 2 + and <u = nw, 

(n = ±1, ±2, ...). This implies that, from comparing the result 
with cold plasma limit, the very large B q plasma behaves just like 
a cold plasma except for the intrinsic cyclotron harmonic behavior. 

Wow in order to get some idea on the relation between the 

equilibrium distribution and instability of the system, we let 

P o p 

at — (uj^ - u^) + i ( 2u) r «^ ) ■ Then separate the dispersion function 
(3) into its real and imaginary parts, we have 


Rep(k,u>)] = P T 


-A 2 


2 

ID + 

p 


£ 

n=l 


y-, v 2 2 

B (k) n u) 
rr ' c 


L 2 2 \ 
P - n (D 
r c I 


[i> 


2 2 2 2 

n (d + P. , 
cl 1 J 


0 


(7) 
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Im[D(k,(u) ] 



£ n 
n=1 [( 


B n (k) 


O 22 

P r - n <" c 


2 + P 2 ' 

l. 


= 0 


( 8 ) 


where P r « - a|, * 2<VV md V k ) = a n W’ 

Now assume that B n (k) ^ 0 for all n £ 1, then Eq. (8) implies 

that either <d =0 or uu . =0 and there is no propagating instability. 
T 1 

If there are more than one particle species in the plasma, the above 

conclusion is still correct so long as each velocity distribution 

satisfies the condition B (k) £ 0 for all n ^ 1. One class of 

distribution functions that satisfies this condition is defined from 

Eq. (5) by Sf Q /&v £ 0 for all v^ > 0. Hence a necessary condition 

for propagating instabilities is that the inequality d f 0 /d v A > 0 

must be satisfied for some range of v . 

Instabilities associated with the class of distributions 

with ftf /fYv £ 0 have been investigated by many, authors (Hamasaki, 
o' 1 

1968a, b; Davidson and Wu, 1970; Tzoar and Yang, 1970; Bornatici and 

Lee, 1970; Lee and Armstrong, 1971; Oaf fey, Thompson and Liu, 1972, 

1973) • Our interest is to investigate the type of distribution with 

?sf /Av > 0 for some range of v . We will take the ring distribution 
w o' 1 x 

which occurs naturally in the earth magnetosphere when high energy 
charged particles, streaming from the sun, are trapped by the earth's 
magnetic field at the bow shock (Kennel and Petschek, 1966). 
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III. INSTABILITY ANALYSIS 


The equilibrium, distribution is taken to be the ring distri- 


bution 


f - 


fi(v - 
2n v 


v ) 
OJ. 


Ol 


/2^ 


2 \ 


exp 


oil 


1 2 v 


oil' 


(9) 


The distribution function is normalized so that 


f f dv = 1 
J o 


( 10 ) 


The dispersion relation can be -written from Eqs. (3) and (6) as 


D(k,u)) = 


cu 


a 


an 2 - p 


Po - P ll p o 


E 

n-~co 


n J 2 (p ) 


f2 - n 


<0, 


ctfl 2 - £ 


2 

+ 


Vo. 

sin Hn 


dr sin Hr sin t J 2p 
o \ r o 


= 0 


cos 


( 11 ) 


roM 
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where 


and 


kv 


Po = 


o± 


O tt) 


^ <i> 

(1 S , 

0) 

c 


2 

V 

a - > 

c 


P., = 




2 

n mv „ 
o oil 


/ 2 ) 


/ v o») ( B >») 


p 


l#»c) l P o m ' r °i / ' g l 

■ (= 2 /4) ' [*>■ 


If we consider an equilibrium distribution of the counter- 
streaming type, say, 


f 

o 


5(v 


V > 
01 


■(v„ - u)‘ 


-( V M + *Y 


^ V 01 2/2tt v 


exp 


oil 


2v‘ 


+ exp 


oil 


2v‘ 


oil 


( 12 ) 
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the dispersion relation will have the same form as Eq. (ll) with 
replaced by pj = [u> 2 /u> 2 ]/[c 2 /(v 2 , + u 2 )]. Distribution functions of 
the form (12) produce no current along the magnetic field. Thus 
the effect of streaming is only to increase the value of P (j . 

In order to determine the curves of marginal instability in 
the parameter space a - B - we study the threshold conditions 
for the onset of complex solutions of the dispersion relation with 
uu. > 0. The necessary condition (Briggs, 1964) for an onset of a 
pair of complex conjugate solutions of D(k,cu) » 0 at ft = ft, P c = P Q 
is that (Q,p o ) be a double root of ft for p Q = p Q , i.e., 



ft=ft 


p o p o 


(13) 


The pinching -pole criterion (Briggs, 1964; Derfler, 1967 ) 
which is necessary for absolute instability is that D(k,u>) = 0 has 
a double root p Q at ft = ft, p Q = p Q , i.e., 

J£. 

5p o 



= 0 


(14) 
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where ft, o satisfy the dispersion relation 
r o 


D(p o , o, a, P |( , p ± ) = o 


( 15 ) 


for some specific values of O., (3 ^ , 6^. 

Note that it is possible to have an instability with zero u> r 
which corresponds to purely growing modes. By setting ft = 0, Eqs. 
(13) — (15) reduce to 


= 0 (identically zero) , 


5ft 


ft=0 


(16) 


r 5D I 

2u> 2 

= c , 

5p oJ 

a t 


z> 

ji 

0 

r 


P„ - p 


p -p 
K o r o 


+ J l ( Po } " J l ( Po )] 


= 0 , (17) 




B(P 0 . 0 - 0) - ^ t-P, - f a + »,P. J o ( y J A )] - 


= 0 


(IB) 


Equations (17) and (18) give a set of two simultaneous equa- 


tions , 
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1 
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VPq> Vp„> 

A/ 

p 0 


(19) 


( 20 ) 


The right-hand side of Eqs. (19) and (20) are plotted in Figure 1. 
Since p o is real, Eq. (19) can be satisfied only if P (| > 1 and for 
a given value of P , Eq. (20) can be satisfied only when the value 
P /P^ is rendered such that the hyperbola [l/P^][l + (P^/p^] = const, 
intersects with the curve [2J q (p q ) J 2_(p 0 )Vp 0 = co 118 ^* When this 
occurs there exist purely imaginary roots of the dispersion relation, 
implying that plasma fluctuations will grow in time without propaga- 
tion. The situation when these two curves just touch gives the 
curve of marginal instability for purely growing modes in the P |( - P 
plane and is presented in Figure 2 parameterized as au 0 = 0. Note that 
the curve is independent of a:. 

Numerical study of the dispersion relation has also been done 
to deter min e the curves of marginal instability for the first three 
propagating frequency bands. Solving Eqs. (15), (l^+), uud (15) 
si mult aneously for various values of Cl, p Q , P^ with a, P^ as para- 
meters, we have obtained the marginal instability boundaries in the 

p - p plane with ot as parameter, as shown in Figures 2 and 5. It 
1 A 

is interesting to note that these instability boundaries are 
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essentially straight lines. Points lying on the right-hand side of 
these straight lines are unstable. Figure 2 is obtained for a = 0.1. 
Note that these three lines, which correspond to the first three 
propagating frequency bands, intersect each other. The effect of 
the perpendicular temperature is destabilizing; increasing the value 
of a decreases the value of P (| for instability, as shown in Figure 5. 
Attention is also called to the fact that the instability boundaries 
corresponding to different values of ct for the same frequency band 
are parallel. Instability can occur for very low 
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IV. D IS PEES ION CHARACTERISTICS 


In order to make sure that these curves truly represent the 

instability boundaries, we have to solve the dispersion relation (ll) 

numerically. The dispersion characteristics for the points A, B, 

C. D, E, F, G, and H in p - P plane for a = 0.1, as shown in Figure 
’ i II 

2, are shown in Figures 4(a)— 4(h), respectively. Cutoff is ob- 
served at each harmonic of the cyclotron frequency and at the plasma 
frequency The mode which starts from is well approximated 

by the cold plasma wave u > 2 = u ) 2 + A 2 . Resonances are found at 

«n harmonics of cu . These observations are consistent with the 

c 

analysis of section II. 

2 2 2 3 

Now we assume u> ~ nu> and (ouu/u)_)(v /v ) « 1, then the 

c pc 0 1| OJ, 

infinite series in Eq. (11) can be well approximated by one term and 
the solution of (ll) is approximately given by 


<u « nu) c 


] 


2 


oil 


OJL 

2 2 

n tu c 


d 

d Pr 


2,2 

- c k 



( 21 ) 


which accounts the undulation of each mode near each harmonic of the 
cyclotron frequency. This is shown in Figure 4(a). It will be noted 
that Eq. (21) predicts that the modes pass through points defined by 
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the zeros of J (p ) and (d/dp)J (p ), in agreement with the exact 
n o n o 

numerical solution. In order to prove that this is true in any case, 
it is only necessary to obtain an alternate representation of the 
dispersion function (ll) (see Appendix B) 






sin Qtj 



( 22 ) 


Clearly if H is an integer, Eq. (22) implies 


J n<"o> 




■ 0 


verifying that the modes pass through harmonics of the cyclotron 
frequency when p Q is a zero of the nth order Bessel function or its 

first derivative. 

o p 2 2 

As (m /m )(v ,,/v ) increases, the loops above a given harmonic 

jy c or oj. 

approach the loops below the harmonic immediately following it. The 
points at which the loops can couple must always lie between 7 R m 
and 7 n+1 m , where m represents the mth zero of J^(p Q ). At point 
B of Figure 2 the n = 1 mode touches the zero frequency axis, thus 
there is a range of p Q in which purely imaginary solutions exist. 

The imaginary solutions are indicated by a broken line, as shewn in 
Figure 4(b). 
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At point C of Figure 2 there are two unstable frequency bands . 
In addition to the purely growing frequency band there is a complex 
frequency band with cu c < tu r < 2ou c . The real parts of the frequencies 
are indicated by fine lines and the corresponding imaginary parts are 
shown by broken lines in Figure 4(c). The presence of these complex 
solutions has the important practical significance that an individual 
propagating mode will grow in time to an amplitude limited only by 
the validity of the linear theory which has been used in obtaining 
the dispersion relation. 

The points D, E, F, G, and H of Figure 2 give the consistent 
results which are shown in Figures 4(d) — 4(h), respectively. It is 
pointed out that the real parts of the complex frequency bands are 
centered at about (n + l/2)a> e and the maximum growth rates of the 
unstable waves decrease with increasing w and are typically of the 
order of the cyclotron frequency. The growth rate can be very strong 
for some appropriate choices of plasma parameters . 

If we further increase the values of f3 „ and (3 , we reach 

li l 

the situation that there are two purely growing bands; they couple 

at o = p ' (a. {3 , p ) and for p > p ' there always exist complex 
r o o II i r o o 

solutions with m r starting from zero at p Q = to infinity as p Q 
approaches to infinity. This situation is typically represented by 
Figures 5 . 
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V. CONCLUSIONS 


The instability of the ordinary electromagnetic mode propagat- 
ing perpendicular to an external magnetic field has been investi- 
gated for a single-species plasma with ring velocity distribution. 


Sf (v (| , v )/jiv > 0 for some v^ > 0 is a necessary condition 
for the existence of propagating instabilities. The instability 
criteria give the marginal instability boundaries not only for the 
purely growing mode, but also for the propagating growing modes. 

These instabilities are all absolute because they satisfy the abso- 
lute instability criterion. The necessary condition for the purely 
growing mode to set in is P^ > 1 and the instability boundary for 
purely growing modes depends on and P^ only, independent of 
a fa = v 2 /c 2 ). The marginal boundaries for propagating instabilities 

\ Oj/ 

axe straight lines in - P plane with P > 1 ^ as a parameter. 

For different a the marginal boundaries of propagating modes with 

u) in the same harmonic branch are parallel. Increasing a decreases 
r 

ft required for instability. For appropriate choice of parameters, 
growth rates may be of the order of the cyclotron frequency. The 
growth rate is enhanced by increasing P . The investigation of loss- 
cone type velocity distribution and the case of two-species plasma 
for the ordinary EM modes has been underway and the results will be 
published in the near future. 
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APPENDIX A 


From Eq. (5), with fl = u/u> c > we have 


„ . ® 2ttcju n r “ o 

D (k,») = » 2 - A 2 - + E / v » dT 

- r \ — — cc > -go 


ns=-oo 


r ” df p 

X / dV j. aV J n (p) 

J o 1 


2 

OD - C 


^ - 4 - 2 ”p / av » v » / ipf < 


nJ„(p) 


Sp 


r- n vr / 
n - n 


n=-° 


(A-l) 


By using the identity £ J^(p) = 1 and ^ he irrtegra3 ' re P re ~ 


n=-c° 


sentation for the Bessel function 


^(p) 




sin 


cos n-r dT 


(A-2) 


the sum of Bessel functions in (A-l) can be written as 
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E 


n=-» 


nJ^(p ) 
2 - n 




n=-«. 


j 2 (p) 

n M-/ 

2 2 
n - n 


1 


1 

n 




cos tit 



1 


(A-3) 


If 2 is complex we have, by means of a Soramerfeld-Watson transforma 
tion, 


cos ht 1 

— a 2 2 = 2i 

n=-«° 0 . - n 


E 


dz 


COS Zt 


c 1 + c 2 t! “ M c 2 


(A-4) 


where the contour + C 2 is taken around the real axis in the 
complex z planes as shown in Figure 6. 
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f 

Jc i^2 


dz 

tan zrr 


cos zt 



£ 


C l +C 2 


dz cos (tt ~ t)z 
sin ZTT n 2 - z 2 


£ 


dz 


sin rz 


2 2 
C ± +€ 2 ft - z 


(A-5) 


The second integral in (A-5) vanishes because its integrand is 
analytic inside the contour C. Since on the large circle 
the contour integral approaches zero as the radius R approaches 
infinity, we have 


r 

o 

Jc 


dz cos (tt - t)z 


C l +C 2 


sin zjt 


ft 2 - 


' C 1*R, 


dz cos (tt - t)z 


sm zrr 


n 2 - 


Jc 2^r, 


dz cos (tt - t)z 


sm Zrr 


n 2 - 


2rri cos (tt - T~)ft 
ft sin rr2 


(A-6) 


Hence (A-4) becomes 
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E *° cos nT _ rr cos (tt - t)z 

J2. 2 _ fi sin fi-rr 

n=-«= fi - n 


(A-T) 


If fi is real we have 



cos nT 



1_ 

2i 


i 


C l +C 2 


dz cos zt 
tan ZTT ^2 _ z 2 


+ 


TT COS 

Cl tan Hit 


(A -8) 


and 


1 

2i 


l 


C l +C 2 


dz cos zt 1_ 

tan z n ^2 _ z 2 2i 


r 

Jc 


C l +C 2 


dz cos (tt - t)z 
sin ^ fi 2 - z 2 


1 

2i 


l 


dz 


C l +€ 2 


sin zt 



. (A-9) 


The first integral in (A-9) vanishes because the integrand has no 
poles in the complex plane except on the real axis . Thus 


/ 

J °1 ¥C 


2 


dz 

tan zrr 


cos ZT 



fi sin aT 


2i 



25 


Hence (A-8) becomes 


c"» cos nr _ tt cos (tt - t)z 

r, « o 2 n 2 ~ a Sin 
n—- 00 u - n 


(A-10) 


which is analytic for all ft except possibly at ft = 0, ±1, ±2, ... 
Note that (A-10) has the same form as (A-7). Hence (A-3) becomes 



nJ 2 (p) 
n K _ 

ft - n 



x \ ft cos (tt - T^ft 
2 I sin nft 


1 . 


(A- 11) 


Finally Eq. (A-l) becomes 


D(k,») - - A 2 - 4 + 4 % / dv » v « / dpf c 


dT J x f 2 P sin 2) 


ft sin o cos (tt - r)ft 

2 (A-12) 

sin ftrr 


Now put rr - t — 9 • 
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j d.T J ]_( 2 P sin sili 2 003 ~ 

-S o 


= J de J x (2p cos |j] a cos ^ cos 20 


= sin 29 cos ^ J^|2p cos 


TT 


0=0 


-If 

J o 


de sin [ 2 P cos § J l( 2p 


cos | 


-if 


d0 sin fte sin 


e a 


2 s(2p cos | 


^2p cos ^ J 2p cos 


£ 

2 


f de sin 20 sin 0 J Q [ 2 P cos |>) • (A-13) 

■Jo 


Thus with the aid of (A-13) the dispersion function reduces to the form 


2 2 2 2 2 f“ 2 f” dp pf c 

D(k,m) - « - c k - w p + 2nco p ^ ▼, J ° iirTo^ 


J dr sin 2t sin t *J q ( 2 p cos -g 


(A-14) 


io|a> 
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APPENDIX B 


From the identity (Watson, 1958) 


2 

TT 



(2p cos t) cos (p, - v ) t dr = J (p) J (p) 

(i- V 


(B-l) 


If we put p, = -v = ft, then 


2 f V 2 

(p) J ^(p) = - J J Q ( 2 p cos T ) cos ft T d T (B- 2 ) 


ft 


and 


#[Vf> 


-2 r 77 / 2 

— / j^(2p cos cos T cos 221 

77 Jo 


(B-3) 


Putting 0 = t/ 2, we have 





cos 75 cos nede 


(b- 4 ) 
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By using Eq. (A-13) we have 

£[ J fl (p> J -n<^] " sS £ d6 sin fl8 sln 6 J pl 2p 008 2 ) ' 

(B-5) 


From Eq. (ll) we finally have 




J_ [j (o') 

p 0 - a- b P tv p; 



(B-6) 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


FIGURE CAPTIONS 

Plot of [2J o (p o ) J 1 (p Q )]/p o and J^(p Q ) - J^(p Q ) versus 

/Vi 

Po- 

Instability boundaries in the P - P plane. u> Q = 0 

j I X 

curve is the instability boundary for purely growing 
modes, -which is independent of a. 1 < uu Q /u) c < 2 , 

2 < < 3, 3 < to q_J (w c < 4 curves represent the 

instability boundaries for the first three propagating 
frequency bands with a = 0.1. Unstable regions lie on 
the right of the curves. 

Instability boundaries in the (3 ^ plane for the 

first propagating frequency band with a as parameter. 
Dispersion characteristics of ordinary electromagnetic 
mode for o: = 0.1. 

(a) P =0.2, (5^ = 1.6 represents point A in Figure 2. 
Round points are defined by (tu = nut) c , 

a j ( p )/dp Q = 0) and triangular points by 
(tu = nc« c , J n (P Q ) = °)- 

(b) p = 0.25, p^ = 1.9 represents point B in Figure 2. 
The broken line is the purely imaginary solution. 

(c) (3^ = 0.25, P n = 2.25; 

(d) = 2 . 0 , p n = 2 . 75 ; 
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Figure 5 




\ 

(e) = 7.0, P n = 4.0; 

(f) p = 13.0, p = 5.o; 

1 II 

(g) p = 0.75, P = 4.0; and. 

(h) P^ = 0.75, P,, = 6.0 

represent points C, D, E, F, G, and H in Figure 2, 
respectively. The real parts of the complex roots 
for frequency are indicated by fine lines and the 
corresponding imaginary parts by broken lines . 
Dispersion characteristics for ordinary electromagnetic 
mode -with o : = 0.1, P^ = 10, P^ = 80. The real parts of 
the complex roots for frequency are indicated by fine 
lines and the corresponding imaginary parts by broken 
lines . 


Figure 6 


Contours in the complex z plane for evaluation of the 
integral in Eq. (A-4) . 
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Figure 4(c) 
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Figure ^(h) 
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